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We map spin ladders with n; legs and couplings J' across all rungs and J(l ± 7) along the legs, 
staggered in both directions, to a sigma model. Setting its 8 — (2m + l)ir (where it is known to 
be gapless), we locate the critical curves in the 7 versus *j plane at each m, and spin 5*. The 
phase diagram is rich and has some surprises: when two gapped chains are suitably coupled, the 
combination becomes gapless. With n;,7 and J'/ J to control, the prospects for experimentally 
observing any one of these equivalent transitions seem bright. We discuss the order parameters and 
the behavior of holes in the RVB description. 

PACS numbers: 75.10.Jm, 75.50.Ee, 75.30. Ds 
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With this paper we contribute to the explosive growth 
in the theoretical and experimental studies of antifer- 
romagnetic spin chains and ladders. Let begin with 
Haldane'sO mapping of the spin-S* Heisenberg chain with 
hamiltonian 



H = jJ2 s ( n ) ■ S(n+ 1) 



(1) 



via spin coherent states on to the nonlinear sigma model 
with euclidean action 
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dxdr 
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Here $ is a three component unit vector, and 9, which 
multiplies i times the integer valued winding number W, 
is 



9 = 2ttS. 



(3) 



There is a hamiltonian derivation of this result due to 
AfHeckcl which we shall allude to later. 

Since 9 enters the path integral via e l8W , it matters 
only mod 2tt and when 9 — 0, n, the path integral is in- 
variant under x — ► —x (parity) under which W — > — W. 
One now argues that integer spin chains have a gap and 
half-integer spin chains do not. First, all integer (half- 
integer) spin chains essentially have 9 = (tt). The sigma 
model with 9 =M is known to have exponentially decay- 
ing correlations^. As for 9 = tt, since the spin-i Bethe 
chain is gapless, so must be the 9 — tt sigma model , 
provided the mapping to the sigma model (most reliable 
for large S) is valid down to spin i. Since so much in 
this paper hinges on the masslessness af the 9 = tt model, 
we point out that Shankar and Readel have shown inde- 
pendently that the sigma model at 9 = tt is massless by 
considering the r-continuum hamiltonian of the lattice 
regulated model. 

It is also accepted that nonstaggered 
legs are gapless only if niS is half-integei 



with n/ 
. This is 



most transparent when the interleg coupling J' is much 
larger than the intraleg coupling J, for we can first solve 
the problem of ni spins along a rung, take the lowest en- 
ergy multiplet in each rung, and then couple them with 
J, thereby getting a single chain, about which every- 
thing is known. One then verifies-that nothing changes 
as J' is lowered. EquivalentlyEErO one can show that 
the topological terms for the chains are additive, giving 
9 = 2irniS. Different modifications of spin ladders can 
also be considered by adding next-to-nearest neighbor 
couplingstj'li 3 ! Coupled spin chains have been studied 
by combining mean^eld theory techniques with exact 
results for one chaint-3. 

The spin systems to be considered here have a very 
important feature: they have staggered weights. Let us 
then begin with a single chain for which 



J(n) = J(l + (-1)" 7 ) 



(4) 



is the coupling between sites n and n + 1. Notice that 
7 — > —7 amounts to sublattice exchange n — > n + 1 and 
that the restriction I7I < 1, keeps the interaction antifer- 
romagnetic. ■— .p. 

Affleck and HaldandiSa showed that in this case 



2 7 rS'(l + 7) 



(5) 



so that when 7 is varied from —1 to +1, 9 passes an 
odd multiple of tt i.e., the system is critical, exactly 2S 
times. It is instructive to interpret these transitions 
the valence bond terminology of Affleck et al. (AKLT)E3, 
wherein each spin-S is viewed as a symmetrized product 
of 25* spinors. As 7 is raised from —1, the chain goes from 
being fully dimerized with all the valence bonds (spinor 
contractions) on odd-n links, to being dimerized with all 
valence bonds on the even-n links. As each spinor index 
switches loyalty, it necessarily reaches a point when it 
can equally well go either way, producing a nonstaggered 
i.e., a gapless spin-^ chain. (It is important to realize 
that the effective interaction of these spin-i degrees of 
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freedom can be nonstaggered even though the original 
Heisenberg interaction is.) 

Consider staggered chains shown in Fig. [l] with hori- 
zontal couplings on the a th -leg (a = 1, . . . , nj) obeying 

J a (n) = J(l + (-l)"+ Q + 1 7 ) (6) 

i.e., staggered in both directions. We now show that such 
systems have a rich phase structure in the 7 versus J' / J 
plane at each ni and S. 

Recall Affleck's derivation of the sigma model hamilto- 
nian from the spin chain by pairing spins, forming their 
difference and sum, and turning these into the sigma 
model field and its conjugate momentum respectively, in 
the li»it of large S. This method was generalized by 
SierracJ to uniform ladders. The main difference was the 
n;-fold increase in the number of degrees of freedom due 
to the transverse label a = 1 , . . . ni for the legs. A low 
energy analysis indicated that only one of these modes 
remained low in energy and defined the effective sigma 
model, while the rest had a gap of order J'. This effective 
model had 9 = 2irniS, (independent of couplings), yield- 
ing the previously quoted result for nonstaggered ladders, 
namely that only an odd number of half-integer chains 
were massless. 

We extended this derivation to the staggered ladders 
and found 

6 = 27rSn l (l+ 1 f ni (J'/J)) (7) 

where 

fm(J'/J) — — T-^rti.odd 
f 

ra=i,3,...,nj-l sm ^2nJ L ^ J cos ^2nJ 

with 5 ni . oc id equal to f if Ui-4s odd and zero otherwise. We 
refer the reader to Sierraliil for a very similar derivation 
in the uniform case. 

The critical points follow from setting this 9 equal to 
odd multiples ofn. 

We analyze two cases: n\ = 2,5 = i,f which should 
convince the reader of the soundness of this method, and 
facilitate the discussion of the cases with larger 5* and n\ . 

Consider Fig. ||. On the 7 axis, where the chains 
decouple, there is just one (which equals 2S) critical point 
corresponding to the nonstaggered spin-i chain. Our 
theory predicts that as we turn on J 7 , this becomes two 
critical points that move towards the walls I7I = f . It also 
tells us that although staggering or interchain coupling 
are individually bad for criticiality, a certain combination 
can sustain criticality. Can we believe this? Let us go 
to 7 = — 1. Now each chain breaks up into disconnected 
pairs, but the disconnected pairs of one chain do not lie 
opposite those of the other, but displaced by one unit. 
When these get coupled by J', we have a "snake" chain 
that winds through the lattice. It is a spin-^ chain with 
staggered weights 2J and J'. Clearly at J' = 2J, it 



becomes critical as predicted by the theory. Thus the 
vertical J' axis is seen to play the role of an effective 7 
for the snake. We display this by showing three snakes on 
the left margin of Fig. |2j, with vertical bonds which are 
stronger than, equal to and weaker than the horizontal 
ones ( 2J). As expected, the same thing happens on 
7 = +1, with n — > n + I . 

Although it is not so easy to understand criticality as 
we go into the rectangle, by continuity of 9, the criti- 
cal curve must exist. There is however one caveat: the 
phase diagram in Fig.pl does not strictly follow from the 
equation for 9 when J —* 0: the two critical curves com- 
ing down from J' / J = 2 on | 7 | = 1 will cut the 7 axis 
at distinct points on either side of the origin instead of 
meeting there. But we know that the sigma model map- 
ping is doomed to fail as J' -> 0: we will get not one 
low energy field, (the putative sigma model field) but 
two, since the gap that separated the sigma model field 
from the other, of order J', vanishes. Fortunately, on 
the 7 axis, where the chains decouple, we know every- 
thing: there is only one transition at 7 = which the 
two chains undergo simultaneously. Thus the final phase 
diagram was obtained by combining what we know on 
the 7 axis (about decoupled chains) with what we know 
off the 7 axis (from the sigma model)tZl. Even off the 7 
axis the sigma model is only to be taken as a guide to 
the topology of the phase diagram and not for the exact 
location of the critical curves. This is because the for- 
mula for 9 is generally not exact except when 7 = and 
9 = 2nniS, in which case the sigma model is invariant 
under parity and you cannot alter 9 by a small amount 
( say of order I / S ), without violating parity. 

Along an arc starting at 7 = — 1, J' = and ending at 
7=1, J' = 0, 9 rises continuosly from to 47r. The criti- 
cal behavior is the same across any of these critical curves 
and the gap will behave as i 2 / 3 , where t is the control pa- 
rameter, as predicted by Cross and FisherES. (There will 
be logarithmic corrections since the 9 = it sigma model 
differs from the conformally invariant WZWpjaiodel by a 
marginally irrelevant operatoru. Chitra et ali3 avoid the 
log by adding a special value of nnn coupling to the spin 
— 2 chain and find an exponent very close to 2/3.) 

Let us examine Fig. [| in tyrms of the RVB picture of 
White, Noack and ScalapinoQ for a nonstaggered n\ = 2, 
S = h system which corresponds to a point vertically 
above the origin in Fig. ^, with some generic J' . (It might 
help to consult their Fig. 3 for this discussion.) In the 
absence of defects, the bonds in each two by two square 
resonate between being vertical (with coupling J') and 
horizontal (with coupling J). A defect forces the bonds to 
be horizontal, staggered, and nonresonating, till we reach 
the next defect. This causes a linear confining potential 
and restricts the excitations to spin-I. In our problem, 
we are free to move to the left of this point towards nega- 
tive 7. Now the staggered horizontal bond configuration 
between the defects becomes more favorable and we soon 
hit the critical curve on which the staggered bonds con- 
figuration becomes degenerate with the resonant ones, 
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and the defects (spinons) are liberated. To the left of 
the critical curve confinement resumes, for reasons best 
understood if we drop vertically from the critical point 
to the J' = axis. Now we have decoupled staggered 
chains. The bonds are dimerized in the preferred sublat- 
tices. A pair of defects now forces singlets on unfavorable 
bonds in the region in between. When J' is turned on, 
between the defect, the bonds can resonate since the de- 
fect has lined them up across each other. Increasing J' 
improves resonance and we finally hit the critical curve. 
In general all our critical curves may be characterized as 
those on which the defects are unconfincd. 

What about the order parameter for the different 
phases? Once again it is best to move up the 7 = — 1 axis, 
where we see that the valence bonds go from being hori- 
zontal to vertical. This is just the Affelck-Haldane trans- 
fer of bonds on a chain, but along the length of the snake, 
wherein even/odd bonds turn into vertical/horizontal 
bonds. 

Consider Fig. || for spin-1. Once again on 7 = —1 we 
get a spin-1 snake, which becomes gapless when its stag- 
gering equals ±1/2 according to the sigma model (0),e9. 
Once again the ratio of couplings J' /2J, determines the 
effective staggering along the snake, not to be confused 
with the original 7 for the ladder. Setting J'/2 J equal to 
(1 ± |)/(1 =F 5) we get critical values J'/2J = 3, 1/3. It 
is clear that we can adapt .the nonlocal order parameter 
of den Nijs and Rommelseta (rendered along the snake) 
to describe the Z2 symmetries. We do not discuss spin-1 
further since it resembles spin-^ in other respects. 

For larger values of spin and ni , each single-chain tran- 
sition on the 7 axis splits into n; transtions as we turn on 
J'. The critical curves bend towards the wall (I7I = 1) 
nearest to them. The parameter 9 rises continuously from 
to 4tttiiS as we follow the arc shown in Figure [| There 
are however some differences. First, we get honeycomb 
ladders instead of snakes for larger n;. Next, we no 
longer have an easy way to see the sigma model is even 
qualitatively correct when it locates critical lines for us. 
However, we expect the model to be weakest when ni or 
S is small. Having passed the test there, it seems im- 
mune to further jeopardy. Finally, if niS is half- integer, 
an odd number of lines will emanate from the origin, one 
of which will go straight up to J' = 00 (corresponding 
to nonstaggered odd-n; half-integer spin chains ladders, 
known to be gapless). 

What does the sigma model have to say about holes? 
It was shown by ShankaiO, that in the large S limit, 
the holes in the single chain may be represented by spin- 
less fermions that couple to the sigma model field via a 
gauge interaction. At finite doping, the fermions render 
the 9 term ineffective, wiping out the sharp distinction 
between integer and half- integer chains. The extension of 
this calculation to ladders will tell us if all ladders have 
exponential decay upon doping. 

To summarize, we have considered the phase diagram 
of ladders with staggered couplings by mapping the lad- 
ders to a sigma model and setting its topological coef- 



ficient 9 to an odd multiple of tt (when the model is 
known to be massless). There were a few surprises: we 
have examples here wherein coupling gapped chains leads 
to gapless chains. This is because there is an interplay 
between staggering and interchain coupling which sepa- 
rately destroy gaplessness, but together can conspire to 
keep the system gapless. Thus two spin-^ chains with 
small staggering and small J' can remain massless. At 
all these phase transitions the gap will vanish as i 2 / 3 , (up 
to logarithms) . We expect gapped states to exhibit linear 
confinement of a pair of defects and possibly pairing, if 
the doping is macroscopic. 

It will be worth confirming these predictions by Monte 
Carlo, Density Matrix Renormalization Group, series ex- 
pansions and so on. The sigma model complements these 
approaches: it does not do so well numerically, but man- 
ages to give at one stroke the phase diagram for any 
choice of S and n/. For instance, we know that on 
7 = ±1, where we have a honeycomb ladder, each transi- 
tion of a single chain gets tranformed into n; transitions 
as J' is varied. This accumulation of critical points (for 
any spin, half-integer or otherwise) facilitates extrapola- 
tion to the ordered state in d = 2, although we cannot 
raise n\ too much. ■— ■ 

We leave it to the ingenuity of the experimentalistsEHl 
to find ladders wherein bonds alternate in both direc- 
tions (^), and either 7 or J'/ J, or both, can be varied 
at least slightly. There is also the option of studying the 
honeycomb ladder, an extreme case of bond alternation 
(7 = ±1). Once any such a ladder is found, it will have 
many transitions, whatever be the spin. For example a 
honeycomb ladder with four legs and spin 1 will have four 
transitions as J' varied, say by applying pressure. 

Given the many formal similarities between the 9 pa- 
rameter of the Quantum Hall Effect and the one we have 
here, one can expect strong parallels between coupled 
chains and coupled Hall planes. 

Note added in proof: The phase diagrams have a nice 
extesnion to J' < 0, for the case where J a (n) — J(l + 
7(— l) n ), i.e., the staggering is only along the leg but not 
along the rung direction: if we lower J' from to —00 
each transition point of the decoupled spin-S chain splits 
into ni lines and all 2niS of them flow down to J = —00 
and terminate at the gammas corresponding to the niS 
transitions of the spin 2S chain. 
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FIG. 1. A typical staggered ladder with staggered cou- 
plings J(l ± 7) along the horizontal legs and J' along the 
vertical rungs. 



FIG. 2. Phase diagram for staggered ladders with ni = 2 
and spin | . The solid lines represent critical lines correspond- 
ing to the ^-parameter being an odd multiple of n, (the dashed 
line is just indicative for counting the number of critical lines 
between 7 = — 1 and 7 = +1.) On the margins we show the 
snake patterns associated to criticality (non-staggered) and 
to gaped phases (staggered). 



FIG. 3. Phase diagram for staggered ladders with n; 
and spin 1. See Fig. 2 for similar explanations. 
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